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12 Abstract: We present simulations of space weathering effects on ice deposits in 

13 regions of permanent shadow on the Moon. These Monte Carlo simulations follow 

14 the effects of space weathering processes on the distribution of the volatiles over 

15 time. The model output constrains the coherence of volatile deposits with depth, 

16 lateral separation, and time. The results suggest that ice sheets become broken and 

17 buried with time. As impacts begin to puncture an initially coherent surficial ice 

18 sheet, small areas with a deficit of ice compared to surrounding areas are formed 

19 first. As time progresses, holes become prevalent and the anomalous regions are 

20 local enhancements of ice concentration in a volume. The 3-D distribution is also 

21 heterogeneous because the ice is buried to varying depths in different locations. 

22 Analysis of the coherence of ice on 10 cm scales predicts that putative ice sheets in 

23 anomalous radar craters are < 100 Myr old. Surface frost becomes homogenized 
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within 20 Myr. The simulations show the data from the LCROSS impact and 

25 surrounding region are consistent with the ice deposit in Cabeus being >1000 Myr 

26 old. For future in situ analysis of cold trap volatiles, a horizontal range of 10 m is 

27 sufficient to acquire surface-based measurements of heterogeneously distributed 

28 ice. These results also support previous analyses that Mercury's cold traps are 

29 young. 

30 

31 1. Introduction 

32 A question of both scientific and exploration interest is how water ice and 

33 other volatiles are distributed in permanently shaded regions (PSR) near the lunar 

34 poles. Practically, knowledge of volatile contents and their distribution in PSRs will 

35 aid in mission planning and design for future human and robotic lunar missions 

36 because those volatiles represent a significant potential resource. Scientifically, the 

37 PSRs potentially hold volatiles deposited over the last 2 billion years (Arnold, 1979; 

38 Siegler et al., 2011). This is an unparalleled, accessible record of volatile flux in the 

39 inner solar system. 

40 Analysis of the possibility of water and other volatiles existing in PSRs has a 

41 long history (e.g., Urey, 1952; Watson et al., 1961; Arnold, 1979). Yet, the first 

42 spectroscopic confirmation of their presence inside PSRs came with the impact of 

43 the Lunar Crater Observation and Sensing Satellite (LCROSS) into Cabeus crater 

44 near the lunar south pole (Colaprete et al., 2010; Gladstone et al., 2010; Killen et al., 

45 2010). Lunar data from the 1990s and 2000s provided evidence consistent with the 

46 presence of ice in the lunar polar regions; yet it was inconclusive. Neutron 
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spectroscopy from Lunar Prospector returned data consistent with the presence of 

48 water ice in the near-subsurface of the Moon in permanently shadowed regions by 

49 measuring heterogeneously distributed hydrogen concentrations at both lunar 

50 poles (e.g., Feldman et al., 1998; Eke et al., 2009). Clementine and ground-based 

51 radar returned tantalizing, but inconclusive evidence of ice in lunar PSRs (Nozette et 

52 al., 1996; Campbell et al. 2006). Later, Mini-SAR and Mini-RF instruments on 

53 Chandrayaan-1 and the Lunar Reconnaissance Orbiter (LRO) detected a radar 

54 signature consistent with water ice in some polar craters on the Moon, but not all 

55 PSRs (Spudis et al., 2010; Neish et al., 2011). Similarly, Lunar Exploration Neutron 

56 Detector (LEND) on LRO also detected a heterogeneous distribution of hydrogen 

57 among lunar PSRs (Mitrofanov et al., 2010). In addition, Lyman Alpha Mapping 

58 Project (LAMP) on LRO detected far ultraviolet (FUV) spectra consistent with of 

59 frost on the surface of some permanently shadowed regions (Gladstone et al., 2012). 

60 Various measurements that are directly and/or indirectly sensitive to ice in 

61 PSRs inherently have different sensitivities, fields of view, and resolution. Data from 

62 different instruments often appear to be in conflict about the distribution of 

63 volatiles they imply. E.g., the weakest spectral feature from LAMP was associated 

64 with Shoemaker crater (Gladstone et al., 2012), which had the strongest hydrogen 

65 feature from LEND (Sanin et al., 2012). The impact of LCROSS into Cabeus released 

66 water and other volatiles (Colaprete et al., 2010; Gladstone et al., 2010; Killen et al., 

67 2010; Schultz et al., 2010), but abundances were higher than the background 

68 amounts detected by neutron spectroscopy, implying heterogeneity within that PSR 

69 either in area or in depth (Mitrofanov et al., 2010; Elphic et al., 2011). 
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To add another level of complication, volatiles emplaced in lunar polar 

71 regions are modified over time by space weathering processes including impact 

72 gardening (e.g., Arnold, 1979). Temperatures below 90 K are too low for thermal 

73 diffusion and sublimation (e.g. Vasavada et al., 1999; Andreas, 2007; Schorghofer 

74 and Taylor, 2007; Paige et al., 2010; Siegler et al., 2011); however, volatiles 

75 deposited in lunar cold traps are subjected to other destructive mechanisms 

76 (Arnold, 1979). It is important to understand the modification processes on several 

77 different scales in order to interpret disparate data in a consistent manner. We 

78 present a simulation and sample the results to provide a self-consistent, multi- 

79 dimensional view of lunar polar volatiles. This work suggests that the lateral 

80 coherence, a property necessary for radar detection of ice, disappears quickly. The 

81 oldest ice deposits, invisible to radar, still can produce a signature observable by 

82 neutron spectroscopy. Spectroscopy in the ultraviolet can provide crucial input for 

83 determining any recent or steady state delivery of volatiles. 

84 

85 2. Model 

86 An ice layer in a cold trap, even if is initially homogeneous, is disrupted by 

87 impact gardening. Impacts can act to remove volatiles as well as preserve volatiles. 

88 As impacts emplace an ejecta blanket, the layer protects volatiles that had been 

89 exposed to the surface from losses to sputtering and sublimation. However, impacts 

90 vaporize material, thereby releasing volatiles at the point of impact. As each 

91 location on the Moon has undergone a unique impact history, this process 

92 introduces heterogeneity into the system. 
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93 We incorporate techniques developed to analyze Apollo drill cores (Arnold, 

94 1975; Borg et al., 1976} to create a Monte Carlo model that examines the evolution 

95 of ice in a permanently shaded region on the Moon (Crider and Vondrak, 2003a; 

96 2003b) and Mercury (Crider and Killen, 2005). Previously, the model simulated the 

97 evolution over time of ice with depth in a single vertical column of regolith. We 

98 have expanded the model to investigate lateral scale lengths of properties of ice. 

99 Instead of simulating a single column, we now simulate two columns with 

100 adjustable lateral separation. We consider impacts and gardening in the area 

101 including the pair of columns and allow any impact to affect whichever simulation 

102 columns are in the crater radius or ejecta blanket. Correlations between columns 

103 provide a look at how features persist in both time and distance. 

104 The model starts with an initial configuration of water ice. Here, we use a 

105 coherent layer of ice on the surface of a given thickness and concentration. The 

106 Monte Carlo simulation implements impacts by randomly generating a set of 

107 discrete impacts using the frequency and size distribution of impacts of m > 1 pg 

108 (Neukum et al., 2001). When a discrete impactor hits, its location relative to each 

109 column of regolith is computed. The model adds or removes material appropriate 

110 for the associated crater. The excavation/ ejecta layer is handled consistently 

111 between the columns. If material is removed from one column and emplaced on the 

112 other, the concentration of volatiles in the excavated material factors into the 

113 concentration of volatiles in the ejecta layer assuming a variable volatile loss from 

114 the impact. The simulations provide a simple approximation of impact cratering 

115 that neglects some important aspects of cratering, including compaction and 
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116 slumping. Smaller impactors are treated as an average accumulation between 

117 larger, discrete impacts (Gault et al. # 1972) because computation time for handling 

118 these discretely increases while the effects of the impacts become smaller. Volatiles 

119 at the surface can be lost during exposure to photolysis or sputtering. Thermal 

120 diffusion is neglected because the rate is small for T < 90 K (Schorghofer and Taylor, 

121 2007). That limits this application to the coldest areas on the Moon. 

122 Each model run presented in this paper starts with an initial ice layer on the 

123 surface that is 10 cm thick and comprised of 100 wt.% water ice. This is more 

124 concentrated than expected, but the results scale with the assumed initial weight 

125 fraction of ice. A set of 1000 pairs of columns was run for each lateral separation (1 

126 cm, 10 cm, 1 m, 10 m, 100 m, 1 km). The water concentration as a function of depth 

127 for each pair is recorded at times 1, 2, 5, 10, 20, 50, 100, 200, 600, and 1,000 Myr 

128 after the ice layer is emplaced. 

129 

130 3. Results 

131 Each Monte Carlo run represents a possible scenario given the inputs to the 

132 program. Running 1000 cases with the same program inputs but a different seed to 

133 the random generator produces a set of profiles that are representative of what can 

134 be expected to occur on the Moon. In general, impacts act to cover the ice layer 

135 because the area of the ejecta blanket is larger than the area of the crater. Thus each 

136 location on the Moon statistically receives more burial events than excavation 

137 events. In any individual run, the ice layer tends to stay intact; buried under a 
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138 thickness of ejecta. The thickness of the ejecta layer varies from case to case; but the 

139 burial rate is 1 mm/Myr on average. 

140 On average, excavation events are deeper than the burial events. 

141 Furthermore, excavation events potentially disrupt the ice layer. Excavation 

142 removes a fraction of the ice deposit. The resulting ice layer is thinner than the 

143 original. The ice concentration in underlying regolith is nearly the same as the 

144 original concentration because of the way the code implements an impact. An 

145 impact that is large enough to penetrate the ice layer will also affect the adjacent 

146 column of material on the lateral scale determined by the depth to diameter ratio of 

147 impact craters. As impacts begin to puncture the ice layer, there will be well- 

148 correlated areas of decreased ice concentration compared to the surroundings. 

149 Both the number of these holes and the size of these holes increase with time as 

150 more and larger impacts become likely. As time progresses, buried ice and holes 

151 become prevalent and the few undisturbed sections of ice are rare. Thus at later 

152 times, the anomalous regions should be expected to have higher ice content than 

153 their surroundings. 

154 Figure 1 shows the heterogeneity as a function of depth and column 

155 separation. The difference between the ice content in each depth bin for column 1 

156 from that in the same depth bin for column 2 is calculated. The standard deviation, 

157 defined here as the sum of the squared differences for each of the 1000 runs divided 

158 by the number of runs, is shown for each depth bin at ages (from left to right) 1 Myr, 

159 10 Myr, 100 Myr, and 1000 Myr as a function of lateral separation. Low standard 

160 deviation means that the ice layer is still uniform (e.g., the dark area at the top in the 
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left plot) or that there is a homogeneous layer that lacks ice (e.g., the dark area at 

162 the bottom of the left plot). The high deviations at the top and bottom of the 

163 enriched ice layer show the depths at which some columns have ice and others do 

164 not. The surface becomes homogeneous with low water concentration (fi Ce < 1 

165 wt.%) at t ~ 20 Myr. Deviations in the ice layer increase with lateral separation in 

166 the 100 Myr case and at later times. Samples separated by >1 m would have very 

167 different vertical profiles for deposits aged > 100 Myr. 

168 The next analysis studies the occurrence of dry regions. In Figure 2, for each 

169 of the 1000 runs of 100 Myr and 1 m lateral separation, a data point appears 

170 showing the average ice content in the top meter of Column 1 vs. Column 2. Four 

171 populations appear. The upper right group includes runs in which both columns 

172 retained a significant amount of ice. The population in the upper left is where 

173 column 2 retained ice, but column 1 lost the ice after an impact. The converse holds 

174 for the constituents in the lower right. The lower left contains the trials where the 

175 ice is removed from both columns. The individual profiles of ice content as a 

176 function of depth are widely varied in the set of runs, reflecting varying burial 

177 depths and thicknesses of retained ice. However, the concentration of the ice, where 

178 it is retained, does not change much from column to column or from the original 

179 value. 

180 The number of points in the upper left, upper right, and lower right 

181 comprises the number of runs in which at least one of the two columns is "wet," 

182 which we define here as having an average of 1 wt. % ice over a given depth. 

183 Dividing the number of points in the shaded area by the total number of runs gives 
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184 the probability of finding at least one wet spot through sampling two columns at 

185 that lateral separation. 

186 These probabilities are given in Figure 3 as a function of lateral separation 

187 and age of the deposit. Note that this simulation assumes that the ice layer began as 

188 a 10-cm thick sheet of pure ice and that thermal diffusion is negligible. At most 

189 ages, the probability is higher with increased lateral separation up to separations of 

190 10 m. Greater separation does not improve the probability. 

191 Next we determine the probability of retaining a solid block of ice (Figure 4.) 

192 The condition chosen for coherence is that there is a large abundance of ice (fi Ce > 

193 0.5) and that the difference between the abundance in the two columns is small 

194 compared to the average value (| fi Ce (A)- fi Ce (B)| < 0.2 < fi Ce >). Coherence on the 

195 size scale on order of the radar wavelength produces a strong radar circular 

196 polarization ratio (CPR) enhancement (e.g., Spudis et al., 2010). For Mini-RF and 

197 Mini-SAR, the relevant scale is ~10 cm. Integrated over the original depths 

198 containing ice, the lateral correlation decreases with age. The probability of 10-cm 

199 size ice blocks remaining as a function of age can be fit by the function Pbiock = 0.904- 

/ 

200 !og(t)/3.108. Detectable ice blocks should remain for ~100 Myr. The ice blocks 

201 become too small by 100 Myr. 

202 

203 4. Discussion 

204 Owing to the spatial resolution of orbital neutron spectroscopy, large area 

205 averages from the model runs are the best representation of neutron data. Neutron 

206 spectroscopy, which is sensitive to about 1 m depth, would therefore be detecting 
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207 ice that has been deposited within the last 1000 Myr. These older deposits have a 

208 high degree of heterogeneity in both depth and lateral distribution. If several ice 

209 layers exist and one is on the surface, the standard dry-over- wet model used to 

210 convert neutron spectra into water abundances might underestimate the water 

211 abundance (Lawrence et al., 2011). Thus, quantifying the ice in a surface layer could 

212 improve models of neutron leakage used in converting neutron spectra to water 

213 abundances. 

214 Surface frost could appear two ways. First, it could be steadily deposited by 

215 volatiles migrating from other areas on the Moon. This case was studied by Crider 

216 and Vondrak (2003a) and is not addressed here. Secondly, volatiles on the surface 

217 can result from a partially buried ice layer, which may have resulted from episodic 

218 sources such as cometary impacts. As space weathering affects the surface more 

219 than it does the subsurface, surface abundances decrease rapidly without additional 

220 input. The model suggests that in the top 1 mm, values are heterogeneous soon 

221 after emplacement as ejecta slowly cover the layer and losses to space occur 

222 through photolysis and sputtering. The heterogeneous period is brief; and the 

223 surface frost becomes homogeneously low within a few 10 Myr after deposition. The 

224 presence of surface frost would likely imply the addition of volatiles at a steady rate. 

225 The concentration of surface frost could be used to constrain the rate of source/loss 

226 to the PSRs. 

227 Lateral coherence disappears quickly on the 20-cm scale according to these 

228 simulations. This implies that ice deposits in anomalous craters in the north polar 

229 region (Spudis et al., 2010) would be young features — less than 100 Myr — if the 
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230 features initially began as a 10 cm thick layer and if thermal diffusion is negligible. 

231 A thicker initial layer would take longer to break apart to the point of losing the 

232 high-CPR radar signature. However, thicker initial deposits require a large source. 

233 The model suggests that ice deposits in the larger cold traps including Cabeus, 

234 which lack a high-CPR radar signature, are older than 100 Myr old or started as 

235 layers thinner than those detectable by radar. 

236 Schultz et al. (2010) concluded that the LCROSS impact excavated a 20-30 m 

237 diameter crater. Timing of the observations of water ice is consistent with the water 

238 being derived from a depth of 2-3 m. Furthermore, Mitrofanov et al. (2010) state 

239 that the impact region is enriched in water by a factor of ~2 compared to the 

240 surrounding area. This study defines “enriched" as both columns containing ice 

241 amounts greater than the average ice amount found from the set of runs. The 

242 likelihood of two columns separated by 10 m being enriched to that degree is 

243 negligible until the deposit is 1000 Myr old (see supplemental material). 

244 Alternatively, the enrichment could occur with depth. Neutron spectroscopy is not 

245 sensitive to ice buried below 1-m of regolith. If the ice deposit in Cabeus is buried 

246 significantly below the depth to which neutrons are sensitive, this modeling would 

247 suggest that the ice deposit is older than 1000 Myr. These findings are more 

248 stringent than the constraint implied by the lack of radar CPR enhancement. 

249 

250 5. Conclusions 

251 Data on many spatial scales indicate that ice in lunar polar regions has a 

252 heterogeneous distribution. Space weathering, especially impact gardening, 
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253 introduces heterogeneity even if an ice deposit is homogeneous initially. This 

254 process takes time, thus the degree of heterogeneity is related to the age of the 

255 deposit. Comparing models to data suggests that the ice deposits in lunar polar 

256 regions are generally old and generally buried. Possible exceptions do exist in the 

257 high CPR crater floors in the northern polar region, which are consistent with 

258 younger ages than deposits in the larger cold traps. Surface frost could indicate that 

259 there are some continually delivered volatiles. 

260 The next steps in lunar polar exploration are to characterize the deposits 

261 with in situ compositional and depth analysis. The heterogeneity due to space 

262 weathering guarantees that lateral variations will exist; and a static lander might 

263 have the misfortune of landing in a dry spot. These simulations suggest that a 

264 lateral range of ~10 m is optimal in the trade between minimizing the probability of 

265 probing a dry spot and minimizing mission duration. These simulations assumed 

266 that the ice layer began as a 10 cm-thick layer. Other initial conditions would 

267 produce different probabilities. Mass spectroscopy would identify constituents, 

268 which likely include multiple sources. With a small number of samples with depth, 

269 one could reliably determine what the peak concentration of the ice deposit was 

270 when it was emplaced, assuming thermal diffusion does not play a role. Although 

271 the age cannot be determined with one or two depth samples, increasing the 

272 number of depth distributions taken can increase the significance of the age inferred 

273 using the burial rate. 

274 The modeling described here has application to other airless bodies, 

275 including Mercury. Mercury's cold traps are inherently different than the Moon's in 
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that radar observations are consistent with every permanently shadowed region 
containing relatively pure ice (Slade et al., 1992). They are interpreted to be young 
(Crider and Killen, 2005). This study supports the conclusion that Mercury's 
deposits are young, assuming that the volatiles have never been thermally 
mobilized. As MESSENGER acquires data, a stronger comparison between the Moon 
and Mercury can be made. 
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362 Figure 1. The standard deviation of the value of the concentration of ice as a function 

363 of depth in one column compared to another column separated by the amount on 

364 the x-axis for an ice layer that was 10 cm thick to start after (from left to right) 1 

365 Myr, 10 Myr, 100 Myr, and 1000 Myr. For reference, the average concentration for 

366 the age is provided to the right of each plot. 
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Figure 2. The integrated ice content in one column versus the other column from a 
set of 1000 runs of 100 Myr and 1 m separation. 
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370 

371 Figure 3. The probability that at least one drill of two drill cores has a measurable 

372 amount of water as a function of the age of the deposit and the lateral separation of 

373 the two cores. For deposits that are younger than 100 Myr, the probability 

374 improves with displacement up to a distance of 10 m. Greater distances do not 

375 significantly increase the probability of finding the ice. 
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376 

377 Figure 4. Integrating the ice abundance over original depth from the model runs 

378 then counting the occurrence rate that the ice abundance is large and is similar to 

379 the ice abundance in the adjacent column provides the probability that an ice block 

380 is retained. For the ice in the top 20 cm, the occurrence rate of 10 cm wide ice 

381 blocks decreases with time. 
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